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Q> ■ is best described by a classical pion field. Such a field might be an isospin 

Q_i' rotated vacuum of the chiral c-model or, as discussed in this work, a solu- 



High energy reactions may produce a state around the collision point that 



tion of the equations of motion resultinng from the coupling of fields of this 
model to quarks produced in the collision. In such configurations all directions 



^ ■ in isospin space are allowed leading to a sizable probability of events with, 



essentially, only charged particles (Centauros) or all neutral particles (anti- 
Centauros). (In more common statistical models of multiparticle production, 
the probability of such events is suppressed exponentially by the total multi- 
plicity.) We find that the isospin violation due to the mass difference of the up 
and down quarks has a significant effect on these distributions and enhances 
the production of events consisting predominantly of neutral particles. 



*Electronic address: anselm@lnpi.spb.su 
^Electronic address: mbander@funth.ps.uci.edu 



1 



1. In recent years several authors have suggested that the celebrated Centauro events 
1^, in which no vr'^'s have been observed versus a large number of charged hadrons, might 
be explained by the production at these high energies of a classical pion field; an interesting 
example is the "disoriented chiral condensate" @,|^. The idea is that such a process, con- 
sidered event-by-event, would correspond to the field being along a given Cartesian isospin 
direction. In events where the isospin is oriented (almost) parallel to the 3-rd axis one would 
expect mainly neutral pions while in events where the isospin lies in the perpendicular plane 
predominantly charged pions would be produced. Let (vti, 7r2, vra) be the three Cartesian 
isotopic amplitudes of the classical pion field. As all the orientations are equivalent, the 
distribution in the amplitude n^, is 

dw ~ dn^, tt'^ = nf + nl + ttI = const . (1) 

The number of neutral pions, no, is proportional to 7r| while the total number of produced 
pions, n = Uo + + ~ tt^. With / = no/n, the fraction of neutral pions, one has 
from ([I|), 

this distribution is normalized to unity. 

Obviously @ predicts many more events with a small number of neutrals than do usual 
statistical mechanisms for pion production. In the latter case one expects dw/df to peak 
at / = 1/3 {n^ = n_ = no = 1/3 as 72 ^ oo)) and to decrease exponentially with n as / 
deviates from this value. The distribution (H) corresponds to the limit n oo and gives for 
the relative number of events with the fraction of neutrals less than / 

For a typical Centauro event / ~ 1/100 and P ~ 10%. This seems to be a reasonable number 
as the five "classic" Centauros represent about 1% of events with appropriate energies [0. 

At the other end of the spectrum, near / = 1, the probability of an event having an 
anomalously large fraction of vr^'s is 
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l-P(/) = l-y7~l(l-/). (4) 

We do not have the square root enhancement exhibited in (|^) and instead we find a hnear 
dependence at the end of the spectrum; however, there still is a finite probability of finding 
events with a large number of vr^'s. It is possible that such "anti-Centauro" events have 
been observed |§] and we shall present a mechanism for enhancing their probability over 
that of §). 

The distribution results from exact isospin symmetry. At the quark level this sym- 
metry is rather strongly violated due to the up-down quark mass difference, 7^ m^. In 
this Letter we shall demonstrate that this mass inequality can enhance the probability of 
anti-Centauros. 

2. A class of solutions for the pion field whose dynamics are governed by a non-linear 
chiral Lagrangian was presented in Ref . ||T| . The results of that work may be understood in 
the following simple way. The Lagrangian is 

£ = Il^r [d^Ud^U^) , (5) 

where = 93 MeV and the unitary matrix U is connected to the pion fields by 

r ■ 



[/ = expP-^ . (6) 



For the particular form 

f/ = exp[^r3^(r,^)] (7) 
the Lagrangian (^) leads to the free equation of motion 

d^e = Q. (8) 
For constant unitary matrices Vl and Vr a generalization of (|^) is 

U = VleM^r,e)VR- (9) 
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this is a general class of solutions which has been studied in Ref. |T|. All other known 
solutions I^J^ are particular cases of (H). 

At large distances from the collision point we require the normal structure of the vacuum, 
i. e. U = 1. Likewise we will pick solutions in which 6{r,t) —>■ a.s r ^ oo. This forces 
Vl = Vr and the solutions (|^) reduce to isotopic rotations of (|^. In other words, takes 
the form 



U = exp[iT ■ n 0{x)] 



(10) 



for some direction n in isotopic spin space. A possible scenario for the production of a 
classical pion field discussed in [^J^ is that inside a certain volume around the collision 
point a state corresponding to a constant (in the volume) 9 is produced. This state is 
degenerate with the normal vacuum (in the limit = 0) but is rotated with respect to it 
in isotopic spin space. In this situation is referred to as "disoriented chiral condensate" . 
It follows from (p!0| ) that any solution of (^ describes chiral dynamics. 

We now introduce interactions of pions with quarks keeping in mind that the pion field is 



the chiral phase of the quark field |]TU[. In the presence of pion fields the quark fields should 
be modified 



qL{x) exp[-r ■ne{x)]qL 



qnix) exp[—T ■ne{x)]qR. 



(11) 



The quark mass terms give rise to the quark-pion interaction Hamiltonian 

_ 1 'i 

H = niuuu + niddd ^exp(-T • n 0) (m_,_ + m_T^) exp(-r ■ n 6)q + h. c. , (12) 

where m± = |(m„ ± m^). For the solution (p!0| ) 
Ti. = q{m+ + m^T^)q — (1 — cos 9)q (m+ + m^n^T ■ n) q + sin Oqi'y^ {m^T ■ n + m^n^) q . 

(13) 
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In the normal vacuum ([T3|) accounts for the pion mass term through the existence of the 
chiral condensate {qq) ^ 0. From (|13|) one sees that = —m^{qq)/f^, n = f^^O. 

The distributions, in the parameters 6 and n, of a classical pion field produced in a high 
energy collision are expected to depend on a production temperature T and have the form 

dwr^ Jll dd{x) exp (^-^ j d^x7i)j dn6{n^ - 1) . (14) 

If the quark density in the collision is not too high {qq) should be set equal to its usual 
vacuum value. Expanding around 6 = (0) becomes 

dw^ JYI de{x) exp (^-^^^ J d'xe"^^ dn6{n' - 1) . (15) 



For T = Tc^ 140 MeV |Tl[ and a volume V ~ 100 fm^ the above is exp [-4 < 6"^ >]; large 



values of 9 will not be excited. However, after the functional 9 integration the distribution 
in isospin directions remains uniform leading immediately to (|2D. 

3. Our critical assumption is that in the high density medium created by such collisions 
the quark density and other bilinears in g, q acquire classical values that may be comparable 
to or larger than the vacuum chiral condensate {qq) ^ —(250 MeV)^. From the explicit 
dependence of (|3]) on we see that isospin rotation symmetry is broken. We consider 
two possibilities: either /(x) = {{qr^q)) 7^ or P (a;) = {{qij^q)) 7^ 0, in addition to S{x) = 
{{qq)) sizable. ((■ ■ ■)) denotes the averaging over quantum fluctuations and we 

allow for a smooth (on the microscopic scale) position dependence. The value of S{x) may 
differ significantly from the vacuum value of {qq). 

We first consider the first case, I{x) 7^ 0; although it has less interesting consequences 
it is simpler to analyze. The functional integration over 9{x) in (|T^) (in the quadratic 
approximation) yields 

dw ~ dn^ . (16) 

J\m^S{x) + ■m^I{x)nl\ 

For the dependence of the above on 113 to be significant it is necessary for the second term 
in the square root to be comparable in magnitude to the first one. This is, however, unlikely 
as their ratio is (even for / = n| = 1) 



m.-I{x) m„ — nid {{uu — dd)) 



(17) 



mj^S{x) ruu + rud {{uu + dd)) ' 

with rriu — nid/mu + ~ —0.3 and the second factor less than unity the 713 dependence 
will be insignificant. We reach the same conclusion if we allow other components of grg to 
acquire some classical value. 

The situation is significantly different if we assume that P{x) has a sizable value. Below, 
we shall return to see whether this is feasible, but first discuss the consequences of this 
assumption. We are now asked to evaluate 



dw ~ n d'^i^) I ^^"^ ^-^P [m+S{x) (1 — cos 9) — m^n^Plx) sin 9] | dnS{n^ — 



(18) 



The exponent has a minimum for a non-zero 9 obtained from tan^^ = m-n^P{x) /mj^S{x). 
The functional integral can be done (again in a quadratic approximation) and, aside from a 
pref actor, yields 

1 



~ exp — j d^x +\Jm\S'^{x) + mi.P'^{x)n^ + m^S{x) 



dn^ 



(19) 



Although we could analyze this result it is simpler to consider the situation where 
|m_P/m+S'| < 1. Keeping only the first term in the expansion of the square root we 
obtain (ignoring, in the case S{x) is positive, terms not depending on n^) 



dw ~ exp 

Remembering that / = rig we find 



1 ml PHx) 



d^x 



\S{x) 



dn. 
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(20) 



dw = N{A) 



df 
2v^' 



(21) 



where 



A^^^ld^x 



3 P\^) 
2Tm+J " '^\S{x)\ ' 



(22) 



and the normalization factor 



N-\A) = C dxe^''\ (23) 







Evidently the change in the distribution is important only if A is large enough. In this case 
the distribution ( pT] ) has a minimum at / = 1/2 A and, contrary to the situation described 
by (0), grows as / approaches 1. For A >> 1 an approximate evaluation of ( p^ ) yields 



dw ~ Ae-(i--^)^^ . (24) 
vJ 

This distribution has a peak at / = 1 and is enhanced near that value by a factor 2A over 
that of @) making anti-Centauros more probable. 

We shall now try to estimate possible values for A. While |5'(x)| presumably coincides 
with the quark density p(a;), P{x) can be represented in the form 

P(x) = iR(T{x) ■ Vp{x) . (25) 

Here cr(x) is some spin density, R is a characteristic linear size of the effective volume (or 
characteristic time before hadronization) and ^ is a constant, probably smaller than one. 
Integrating (^) we get 

/ ^^^74 = = o^R —^ P ■ 26 

J p 3 r 

We use r as a characteristic length for the gradient; this variation in density is likely to 
be confined to the surface of the quark matter produced in the collision. We assume that 
the volume over which P{x) does not vanish is AirR^r. The spin densisities are averaged 
approximately to unity. Thus for the parameter A we have: 

^ = ^ ——N ~ T^-^N . 27 
21 m+ r 70 r 

Here = AtiR^p/3 is the number of quarks produced. We believe one could expect N > 200 
in a sphere of -R ~ 3 fm (note that for the vacuum p =< qq >= 2 fm~^, so that ^ 200). 
For R/r ^ 5 we find A ~ 15^^ and for ^ > 0.25 A will be sizable enough to enhance the 
probability of anti-Centauros. Note that ^ < 0.8 is required for the approximation in going 



7 



from ([I9D to (^). We are well aware of the crudeness of these estimates and the purpose of 
this exercise was only to show that values of A> 1 are not excluded. 

The whole change in the distribution of neutrals is due to the violation of isotopic spin 
invariance; the parameter A in ( pT]) is proportional to {rriu — rnd)"^- Can we claim that the 
anti-Centauro events are caused by the mass difference of light quarks? 
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warm hospitality. This research was supported in part by the National Science Foundation 
under Grant PHY-9208386. 
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